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Abstract 

Tliis paper presents new investigations on coherent scattering in the forward direc- 
tion (orders of magnitude; conservation of energy; dependence of scattered Hght on 
geometry and wavelength), and on how scattered light contamination in the spec- 
trum of reddened stars is possibly related to as yet unexplained observations (the 
diminution of the 2200 A bump when the obscuring material is close to the star, the 
difference between Hipparcos and photometric distances). This paper then goes on 
to discuss the fit of the extinction curve, a possible role of extinction by the gas in 
the far-UV, and the reasons of the inadequacy of the Fitzpatrick &; Massa (1990) 
fit. 
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1 Introduction 



When, in the past, observations were limited to the visible spectrum, interstel- 
lar extinction could be understood by the presence of large (compared to the 
wavelength) forward scattering grains, with scattering and extinction cross- 
sections proportional to 1/A; the visible extinction curve {Ax as a function of 
1/A) is linear in 1/A. The column density of large grains between a star and 
the observer is proportional to the slope of the visible extinction A\, which is 
twice the visible reddening E[B — V) = As — Ay of the star if A is expressed 
in /im (Zagury, 2000b). 
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The problem of interstellar extinction was complicated by the release, in the 
1960-70's, of the first UV spectra of reddened stars. Extinction in the UV is 
less than the UV extension of the linear (in 1/A) rise of the visible extinction, 
which can be explained either by a specific composition and size distribution 
of interstellar dust, or by an addition of starlight scattered and re-injected 
into the beam of the telescope. In the former case we do see what we think 
we see, i.e. the direct light from the star extinguished by interstellar dust. In 
the latter the effect of the additional component of scattered light diminishes 
the impression of extinction of the direct starlight by large grains. 

The UV extinction curve exhibits a marked extinction feature, centered close 
to 2200 A, in its near-UV part, and consists of a steep rise, though still under 
the UV extension of the visible extinction, in the far-UV. 

If we observe the direct and extinguished light from reddened stars, the re- 
lation between these features and the large grains, and their variations from 
cloud to cloud, may be investigated by comparing extinction curves in differ- 
ent directions, normalized to a common slope in the visible (i.e. to the same 
column density of large grains) . If extinction curves in different directions are 
normalized by E{B — V), their visible parts superimpose well, but large and 
unpredictible variations in the bump region and in the far-UV are observed 
(Figure 6 in Bless & Savage (1972)): for an equal column density of large 
grains, the observed UV extinction can noticeably differ from one interstellar 
cloud to another. 

Consequently, if the light we receive from a reddened star is dominated by 
the direct extinguished light from this star, large grains do not contribute to 
the UV extinction features, and their UV extinction must be a fiat underlying 
continuum. Different types of particles -a minimum of three- are then needed 
to explain the visible to UV extinction curve (Desert et al., 1990; Greenberg, 
2000). Large grains are responsible for the extinction of starlight in the visible 
and have a grey extinction in the UV. The 2200 A bump feature is attributed 
to a distribution of small grains, and the far-UV rise of the extinction should be 
due to the presence in interstellar clouds of Polycyclic Aromatic Hydrocarbons 
(PAHs: AUamandola et al. (1985)) molecules, which should re-emitt in the 
near-infrared what they absorb in the far-UV. 

The respective proportion of each of these three types of particles determines 
the importance of the extinction in each of the three parts (visible, 2200 A 
bump, and far-UV regions) of the extinction curve. The observed extinction 
curve in one direction is then the sum of these 3 extinctions, and requires at 
least 3 parameters to be fitted. In fact, Fitzpatrick & Massa (1986, 1988, 1990) 
have shown that such a decomposition needs 6 independent parameters. 

However, this standard interpretation of the extinction curve meets observa- 
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tional and analytical difficulties. 



The observed relationships between the different parts of the extinction curve 
(see Figure 4 in Savage et al. (1985), Figure 9 in Savage (1975), CardeUi et al. 
(1989), and Zagury (2001d)), are difficult to reconcile with a three component 
grain model in which the proportion of each type of particles varies with 
the line of sight, and in an unpredictable way. In addition, no dependence 
of the different types of particles on environment has ever been evidenced 
(Jenniskens & Greenberg, 1993), so that the reason for these variations remains 
a mystery. For example, the claim that large grains might be destroyed into 
small grains in low column density, UV-exposed, regions, is contradicted by 
the linear relation which exists (Savage, 1975) between the bump height and 
E{B—V); this relationship means that the column density of small grains must 
grow in proportion to that of large grains. More specific problems remain: these 
include the precise identification of the PAH, the means by which they absorb 
the UV and rc-cmit in the near infrared, and the difficulty of reproducing in 
laboratory particles with such particular properties. 

Recent analysis of UV observations by the International Ultraviolet Explorer 
(lUE) sateUite contradict the near-UV breakdown of the extinction by large 

grains that is assumed by the standard theory. The scattering cross-section 
deduced from the UV spectra of nebulae varies as 1/A (Zagury, 2000a), as 
in the visible, which proves the continuity of the scattering cross section of 
interstellar matter, and a fortiori the continuity from the visible to the UV of 
the extinction law by large grains (Zagury, 2002a). 

The visible- UV continuity of the extinction by large grains was further demon- 
strated for directions where the extinction is very low; the light we receive from 
stars with E{B — V) « 0.1 mag. is extinguished by the same exponential of 
1/A in the visible and in the UV (Zagury, 2001c). At intermediate reddenings 
(~ 0.1 < E{B — V) <~ 0.4) the visible extinction law still extends in the 
near-UV, down to the bump region (Zagury, 2000b). 

Analytically, the situation of the standard theory is no better. Not only have 
Fitzpatrick & Massa (1988) been unable to give any physical meaning to their 
mathematical expression of the far-UV extinction curve, but the Fitzpatrick & 
Massa (1990) fit for extinction curves is curiously confined to the UV domain: 
its extension to the visible systematically diverges from the observations (see 
Figure 4 of this paper and section 12.1 for a detailed discussion on this point). 
Second (and it is stricking that this was passed over in silence in nearly all 
the articles pubhshed so far on the three component grain model) CardeUi et 
al. (1989) have shown that extinction curves can generally be fitted using one 
parameter only, over all the visible-UV spectrum. If one parameter is enough 
to reproduce a large set of extinction curves, why does the standard theory of 
extinction needs 6 parameters? And why, with so many degrees of freedom, is 
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it not able to do as well as the Cardelli et al. (1989) fit, i.e. to reproduce the 
complete visible+UV extinction curve? 

Unless one has absolute faith (Li & Greenberg (2002), 'In Dust We Trust ') 

in the three component model, the several problems met by the standard 
theory of interstellar extinction prompt out to seek for other solutions. The 
only alternative to the standard theory is if the UV light we receive from 
reddened stars is not direct starlight alone, as it has been accepted so far, but 
that it also contains a non-negligible proportion of scattered hght. Scattered 
light will first modify the UV part of the spectrum because it is in the UV that 
extinction, thus the number of photons available for scattering, is the highest. 

The scattered light component is most easily observed in directions of moder- 
ate reddening where it is hmited to the far-UV wavelengths. In these regions, 
it clearly appears as an excess superimposed on the tail of the exponential 
extinction by the large grains (Zagury, 2000b). At large optical depths, direct 
starlight, extinguished by large grains, becomes much smaller than scattered 
light. This is especially true at UV wavelengths, where, if the linear extinction 
by large grains extends to the UV, a moderate column density of interstellar 
matter will be able to extinguish all the direct starlight: we no longer see the 
star but its light scattered by foreground interstellar matter. If the visible 
reddening E{B — V) is important, the scattered light component can be large 
enough to modify the near-UV extinction, and even be detected in the visible 
(Zagury, 2001a). 

The lower level of extinction in the UV, compared to what is expected by Mie 
theory for 'normal' large grains, is due not to a change of the extinction law 
of these grains between the visible and the UV, and to additional particles, 
but to this addition of scattered hght. 

This not only conforms to observations, but also leads to a simplification in 
our understanding of the interstellar medium: large grains follow a continuous, 
linear, extinction law from the near infrared to the far-UV; their properties 
are the same in all directions. The constraints on the nature of large grains 
implied by the standard theory, as well as the need for a grey extinction in 
the UV, are removed. 

My papers published over the past five years have highlighted some properties 
of the scattered light component. To explain both its importance, a probable 
dependence on A as l/A"^, and the small angle within which it is effective, it 
seems that scattering must be coherent, thus confined to a very small angle (of 
order 10~^") around the direction of the star, and due to spherically symetric 
particles (Zagury, 2002c). Hydrogen fuUfills all the conditions imposed to the 
scatterers (Zagury, 2002c). 

In the first part of this paper I will reinvestigate some aspects of coherent 
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Fig. 1. Representation of the angles employed in the text 

scattering from the near forward directions (section 2 to section 7). Amongst 
the consequences, I shaU then discuss the extinction curve of stars with cir- 
cumstellar dust (section 8 to section 10), the spurious effect a scattered hght 
component could have on distance estimates established by means of photom- 
etry (section 11), and the fit of the extinction curve (section 12). 



2 Incoherent scattering 

Let us consider a star of luminosity observed behind an interstellar cloud. 
The cloud is at a mean distance rfg from the star and /q from the observer. 
D = do + Iq is the distance star-observer, d and / are the distance of a point 
M of the cloud, its distance to the star-observer axis is h (see Figure 1). 
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The star is assimilated to a point source of unreddened flux measured by an 
observer on earth: Fq = L^/ {AttD"^) . 



The power per unit surface received by the observer due to the scattering of 
starlight by one atom, with scattering cross section a, at M is: 
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[for sake of simplicity, since we will be interested by orders of magnitude 
only, the Rayleigh cross-section is assumed to be isotropic; this is true within 
a factor of two, see van de Hulst (1969), section 6.12]. 

For an average column density of hydrogen Nh of the cloud, the power scat- 
tered by the gaz reaching earth per unit surface can be estimated by: 

if the scattering is not coherent. 

Considerations on energy conservation show that the power scattered by an 
homogenous medium and received per unit surface will not exceed Fq. This is 
demonstrated by imagining a spherical sheet of gas centered on the star and 
dense enough to extinguish all of its light. The scattered light received by a 
unit surface of a second sphere of radius D will then be at its maximum and, 
for reasons of symetry, P^s = L*/(47rL'^) = Fq. 

Pus, for a cloud half way between the star and the observer {Iq — do — D/2), 
will be: 



^ Stt (/o cos t^~^ j 



L °° 1 



(i + xY 

= 0.5aNHFo (3) 
= 0.5^^ (4) 

with (To = (jA^, and x — tg^O. 

The Rayleigh cross-section of a particle small compared to the wavelength is 
(8/3)7rA;''Q;^ (van de Hulst, 1969), with k — '2,t^/\ and a the polarisability of 
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the particle. For atomic hydrogen, a = {9/2)al = 6.710~^^cin^ (CGS units, 
Oq = 0.52 A is the Bohr radius of hydrogen). The Rayleigh cross-section of 
hydrogen is: a = 5.86 10~^^/A"^ cm^ (A in cm). 

Typical reddenings E(B — y) of a few 0.1 mag. {Ay ~ a few 0.3 mag.) corre- 
spond to column densities of hydrogen Nh ~ 10^^ cm^^ (Bohlin et al., 1978). 
In the UV, 1/A^ is in the range 10^^ — 10^°cm~^. The maximum contribution 
of incoherent scattering by the gaz to the UV light received from the direction 
of the star, estimated from Equation 4, will remain small, 3 10~^Fo at most. 



3 Coherent scattering in the neetr forwetrd directions 

It is only in the near forward directions that scattering is singular enough to 
require a specific attention since the waves scattered by all atoms on the line of 
sight, regardless of their spatial organisation, keep a phase relation (Bohren & 
Huffman, 1983). 

This coherent scattering is superimposed on the incoherent scattering esti- 
mated above, and restricted to a very small vohime of the cloud around the 
direction of the star. As a first approximation, it will be assumed that the 
coherent scattered light can be described by a mean surface 5*0 of the cloud 
around the direction of the star, and by ric atoms of the cloud within 5*0, so 
that the intensity of the coherent scattering the observer measures is simply 
nlPu. 

This assumption separates a focused beam of scattered light (received from 
the near-forward directions), from the classical incoherent scattered light (re- 
ceived from the rest of the cloud). Reality is certainly more complex, with 
successive regions in phase and in phase opposition, as we move away from 
the direction of the star. In a perfectly homogenous medium the contributions 
of these successive zones must cancel to leave only the incoherent component 
of scattered light. The simplified approach adopted here will nevertheless be 
used to derive first orders of magnitude for 5*0, ric, and an upper value of the 
power Pes of the coherent scattered light received by the observer. 

The order of magnitude of 5*0 should be given by the size of the first Frcsncl 
zones, for the wavelength range over which coherent scattering is observed. 
The surface of the first Presnel zone, as it is viewed by the observer, for a 
given wavelength A, is: 

= .A^. (5) 
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Sx is maximum when do — lo — D/2. 



S\ defines the region of the cloud within which all the path-lengths from the 
star to the observer differ by less than A/ 2. It is also the largest surface within 
which the scattered light, at wavelength A, observed at the positon of the 
observer, is constructively coherent. Consequently, if Aj is a lower limit of the 
wavelength range over which scattering is constructively coherent: < Sx^. 

For an atom in Sq, A, the difference between the star-atom-observer path- 
length and D, is (see Figure 1): 

A = = i,^^ (6) 

2 Iquq 2 do 



Viewed from the position of the observer, Sq is seen under an angle 6o = 
{ro/lo)'^-^ {tq is the radius of ^o), which must be less than the angle 6xr. 



Dlo J v20ooAy \ lo J \D 

{9\^ was obtained by setting A = Ai/2 in Equation 6). 

Tic can logically be estimated by NhSq. 

The power scattered by 5*0, at wavelength A, is: 

An upper value of Pes is determined by setting Sq — Sx. 

^cs,max 4A^ 

Note that Pes depends on A~^, while Pcs,max varies as A~^. 
The power Ps^ extinguished by the gas in 5*0 is: 



= aNn-^So = aNnSo^ Fq (10) 
47rao do 



The ratio, Psq/ (rJPcs), of the power extinguished in 5*0 to the power received 
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by a detector of size at the position of the observer is: 



4 Numerical orders of magnitude 



4-1 Numerical application 



To estimate orders of magnitude for the quantities defined in section 3, we 
will suppose do = /q = D/2, and a cloud, of average column density — 
lO^^cm"^ [E{B — V") ~ 0.3 mag.], a hundred parsec far from the sun. The 
raylcigh cross-section of hydrogen at A = 1500 A (1/A = 6/im^^), the UV 
wavelength used in this section, is cr = 10"^^ cm^. The low column density 
approximation applies since aNfj ~ 10~^ << 1. 

The following orders of magnitude are then obtained: 




ric - NhSx ~ 5 lO'^^'atoms (15) 
Pcs,max = ^N^X'^aFo - 2 lO^Fo (16) 

The radius of the first Presnel zone, for astronomical distances, is extremely 
large, of order 500 km, at UV wavelengths and for a cloud 100 pc away. The 
power per unit surface of the light scattered from the first Fresnel zone, at the 
position of the observer, can theoretically be much larger than the unreddened 
flux of the star. 



4.2 Comparison with observations 



If extinction by large grains affects scattered light in the same proportion as 
direct starlight, and if the Rayleigh (exponential in l/A"^) extinction by the 
gas remains small. Pes/ -^0 is equal to the observed ratio of scattered to direct 
starlights. 
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Fig. 2. Ratio, Pcs/Fq, of scattered to direct starlight for a sample of stars. The ratios 
correspond to the j3E^^/X'^ term of Equation 19. They are deduced from Table A.l, 
from Zagury (2001a) for HD46223. 
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Figure 2 plots the ratio of scattered to direct starlights for a number of stars 
for which these two components have been estimated (Zagury (2001a) and the 
Appendix of this paper). The ratio of scattered to direct starlights increases 
with column density, as expected. The effect of gas extinction (the exponential 
of in Table A.l) on the fits of the reduced spectra is small under 1/A ~ 
6//m~^, meaning that for this sample of stars Pcs/Fq is in the range [0.1, 100] 
aroimd 1500 /im. For Nh = lO^^cm^^ [Euv ~ 0.2) the observed Pcs/Fo is 
between 0.1 and 10 (Figure 2), thus substantially smaller than in the idealized 
conditions of sections 3 and 4.1. Pes/ Fq can grow with larger values of E{B — 
V), but it cannot grow too much because the gas extinction term e"^/^ of 
equation 19 will then become important. 



5 Energy conservation 



As long as Fresnel's theory is correct, there is no doubt that NhS\ atoms 
locked in the first Fresnel zone will give a scattered flux, at the position of the 
observer, given by equation 16, and estimated in section 4.1, very much greater 
than the unreddened flux of the source. Because of energy conservation, we 
must have: Pso/{rlPcs) ^ 1- Pso: for a cloud half way between the source and 
the observer, and 100 pc away from the sun, can be estimated by: 4(tNhS\Fo ~ 
4 10"Fo (equation 10), if E{B -V) = 0.2 and for A = 1500 A. From Figure 2, 
Pes is between Fq and IOFq at A = 1500 A (1/A = 6 /xm"^), for E{B-V) < 0.3. 
The condition imposed by energy conservation is then: <^ 2 10^ cm. 

The power per unit surface scattered by NhSx atoms in the first Fresnel zone 
is sufficient to explain the deviation from the linear extinction by large grains 
observed in the UV. However, interstellar clouds are certainly not limited to a 
small space between the star and the observer, and for an homogenous medium 
a unit surface at the position of the observer should not receive more power 
than 10~^Fo (for Nh ~ 10^^ cm"^), and in no way more than Fq (section 2). 
The implication of this apparent contradiction is that interstellar clouds are 
not homogenous at scales larger than Tq , which means that density fluctuations 
must exist at scales smaller than a few hundred kilometers. 

Reality is then in between these two extreme approximations, of a perfectly 
homogenous medium on the one hand, of a concentration of gas in the first 
Fresnel zone (and nothing around) on the other. This is further justified by the 
observational orders of magnitude found in section 4.2 which are significantly 
smaller than the theoretical predictions of section 4.1, and much higher than 
what the brightness of an homogenous medium can be presumed to be. 
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6 Wavelength dependence of the scattered light 

The maximum flux observable from coherent scattering from the flrst Presnel 
zone at wavelength A varies as (Equation 9). However, coherent scattering 
occurs from well defined interstellar structures which each will give a 
dependence (Equation 8), as it is indicated by observation (Zagury, 2001a,d). 

It is also worth noting that the depth of the cloud along the hue of sight 
(which was not considered in the preceding sections) will not change the 1/A^ 
wavelength dependence. As far as (Iq (and Iq) docs not vary appreciably across 
the cloud, the waves scattered by atoms close to the line of sight, regardless of 
their exact projected position on the star-observer axis, reach the observer with 
identical amplitudes and phases, and add positively. If do varies appreciably 
across the cloud the phases at the position of the observer will remain the 
same, though a sum over different amplitudes will have to be considered; but 
the 1/A*^ dependence of the scattered light will be unchanged. 



7 Dependence of the scattered Ught on geometry 

The scattered coherent light from the cloud depends on geometry through 
the dependence of Sq on distances, and through the small scale organisation 
of interstellar matter. For instance, if the scattering medium is close to the 
star, path-lengths' differences will be large even very close to the direction of 
the star; coherent scattering is restricted to a quasi-null surface around the 
complete forward direction (5*0, asS^ —>■ if do —>■ 0, equation 5), and Pes 
should tend to 0. 



8 St£irs with circumstelleir dust 

We expect (section 7) stars with circumstellar dust to have less scattered 
starlight contamination than if they were farther behind a cloud. If, as I believe 
(Zagury, 2002b). the bump is associated with the scattered light component, 
it follows that stars surrounded by interstellar matter should have a smaller 
bump than expected from their reddening. This is exactly the conclusion Sitko 
et al. (1981) reach from a sample of hot stars with a dust shell. 

In the next sections 1 discuss the case of planetary nebulae (P.N.), and the 
shell stars of Sitko et al. (1981) article. Two methods can be used to detect an 
absence of scattered light in the UV spectrum of these object. One way is to 
prove [as it was done in Zagury (2001c) for directions of very low reddening] 
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that the extinction law remains the same in the visible and in the UV: this is 
the case if the UV reduced spectrum of the object exponentially decreases as 
g-2£;„„/A (-j)^ ^j^^ Zagury (2000b)), with E^v ~ E{B - V). Another way 
is to check the presence of a 2200 A bump: the study given in Zagury (2001c), 
on the extinction in very light interstellar media, shows that the threshold of 
E[B — V) above which a bump generally appears should be around 0.05 mag.; 
it is also above this value that the far-UV part of the extinction curve will 
in general no more be fitted by a straight line. The absence of a bump for 
a reddening much larger than 0.05 mag. indicates that the scattered light 
component is negligible. 
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Table 1 

lUE planetary nebulae with no apparent bump 



lUE planetary nebulae with no apparent bump 



name 


\ / 


E(B - V)J^') 


ph 


4t^ 

tv 


lUE file 


NGC6905, HD193949 


0.651 


0.166 


5.4 




LWR10067 


NGC5189, HD117622 


0.6 


b=-3.45 


1.5 




LWR07171 


BD-22°3467, A35 


0.5a 


0.112 




0.134 


LWP22520 


NGC6818, HD186282 


0.217 


0.162 


2.5 




LWR01608 


NGC6720, HD175353 


0.2 


0.095 


400 




LWP20393 


NGC3211, HD89516 


0.20 


b=-4.87 


3.4 




LWR03203 


FB138 


0.2a 


0.074 




0.243 


LWR12596 


NGC2392, HD59088 


0.15 (0.17a) 


0.051 


3.5 


1.1 


LWP24370 


NGC650, HD10346, 


0.14 


0.25 


0.5 




LWP20915 


NGC6853, M27 


0.126 


b=-3.70 


0.4 




LWR05515 


NGC7662, HD220733 


0.105 


0.139 


1.6 




LWP20683 


NGC2371 


0.1 


0.046 


5.3 




LWR04210 


TP94/LS TTT17SQQ1 


U.Uo 


U.lUl 


A A 




T WRn97p;fi 

IjVV IyUZ ( 00 


IC4593, A27 


0.04 


0.058 


3.9 




LWR05705 


NGC4361, HD107969 


0.03 


0.041 


1.6 




LWR15879 


NGC6210, HD151121 


0.03 


0.054 


2.5 




LWR09422 


NGC246 


0.02a 


0.027 




0.63H 


LWR06806 


NGC6826, HD186924 


0.02 


0.11 


1.2 




LWR12580 


NGC7293, Helix 


0.02 


0.292 


0.2 




LWR08942 



(1) Reddening of the P.N.. E{B - l^) ~ 0.7c is calculated from the 'c' (the log- 
arithmic extinction in H/3) column of Table 1 in Cazetta &; Maciel (2000); a 
after the magnitude, indicates that E{B — V) comes from Acker et al. (1998). 

(2) Large scale value of E{B - V) from Schlegel et al. (1998) 

(http://nedwww.ipac.caltech.edu/forms/calculator.html). Gives the red- 
dening due to interstellar cirrus in the direction of the P.N., for \b\ > 5. When 
|6| < 5, the galactic latitude replaces E[B — V)ci- If the P.N. is behind the 
cirrus, E{B — V)ci is the reddening due to foreground dust, excluding the 
material surrounding the star. 

(3) Photometric distance (kpc) from Cazetta Sz Maciel (2000). 

(4) Trigonometric distance (kpc), from Hipparcos, when available, 'oo' is for neg- 
ative Hipparcos parallaxes. 

(5) lUE number of one lUE long-wavel^gth spectrum that can be used to prove 
the absence or the presence of a 2200 A bump. 



Table 2 

lUE planetary nebulae with a bump. See Table 1 for the meaning of each column. 



lUE planetary nebulae with a bump 


name 


HiyjD — V ) 


TT't R T/^ 
-C'l-D — 1/ jci 


0-ph 


atr 


TTTT71 

lurLi me 


xienz-yy 


U. ( o4C 


U.oU 


A 
4.Z 




LiWr^iOUDo 


iNLiUD/Oi, xlUiMDOD 


U. 10 


U.Oi / 


o.y 




T WTT> 1 n77r^ 
LiWrLiU/ (0 


iL'4yy(, xiuiyoooo 


U.IO 


n 1 QQ 
U.ioy 


1 O 




T wm 1 1 m 1 


iNLxvyDD/y, xlUiDy4DU 


U.Do 


U.DUD 


O 1 

z.i 




LiWr^ioozy 


iN*jUDOD(, rUJiDOyoO 


U.oz 


D — -U.DO 


i 




-LW r ZO (Uo 


iNVjL^4U, llUoZD 


u.D (^u.oaj 


n QQ 
U.oy 


i.O 


oo 


-L/VV rv,iOiU4 


JlJJiDMDZ, oWotl 


u.oua 


U.Z ( 4 




U.iiZ 


J_iVV xlU4oU4 


Mnmcii/i wnos2i f;7Q 

iNLT^^1014, xlUzoiD/y 


u.4oa 


U.D ( y 




U.iOD 


T^SARni 070 
J-iVVrtUlZ (y 
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9 Stairs with circumstelleir dust; planetairy nebulae 

For this section I wiU use the large dataset of planetary nebulae (P-N.) from 
Table 1 in Cazetta & Maciel (2000), and the P.N., observed by Hipparcos, of 
Acker et al. (1998). This sample of P.N. covers a large part of class 70-71 lUE 
objects (planetary nebulae observed with and without the central star). 

In Table 1 and Table 2, these P.N. have been classified according to whether 
they have a bump at 2200 A or not. P.N. are presented in decreasing order of 
reddening. 

The striking difference of P.N. with no bump and P.N. with a bump, is that 
the reddening of the former is on the average much lower than that of the 
latter. This increase of reddening of P.N. with a bump compared to that of 
P.N. with no bump correlates with an increase of the 'large scale' reddening, 
Eci, due to interstellar cirrus on the line of sight. Consequently, the difference 
of reddening between the two classes of P.N. (from Table 1 and Table 2) will be 
attributed to a difference of foreground extinction. If this part of the reddening 
could be identified and substracted, it is probable that Table 2 and Table 1 
would be alike. 

Thus, if it is assumed that the two sets of P.N. , P.N. with a bump and fore- 
ground dust from Table 2, and P.N. without an apparent bump and negligible 
foreground dust from Table 1, differ because of the foreground extinction but 
have the same intrinsic properties on the average, it will be concluded that: 

• The reddening of a P.N. due to circumstellar dust is generally low. For most 
P.N., intrinsic E{B — V) should not exceed 0.1-0.2 mag.. Except for a few 
cases, larger extinctions are due to foreground dust. 

• A 2200 A bump in a P.N. is due mainly to foreground dust, not to the 
surrounding material. 

That the intrinsic E{B — V) should be low, 0.1-0.2 mag at most, reflects a 
general feeling, already expressed in the literature (Phillips, 1984). According 
to Phillips (1998), exceptions are compact nebulae of small radius (< 0.1 pc). 
This may be the case of A35, NGC6905 and NGC5189 in Table 1, and of 
IC4997 and HD167362 in Table 2, for which the intrinsic reddening is appar- 
ently large, since the large scale reddening is low in these directions (though 
local enhancements of the foreground column density can not be excluded). 
A35, NGC6905 are cited by PhiUips (1984) (see his Figure 3) as two among 
the few P.N. with higher intrinsic reddening than the average. 

The intrinsic reddening of A35, NGC6905 and NGC5189 should be larger than 
0.5 mag. The absence of an apparent bump in these nebulae, as for other P.N. 
with intrinsic E{B — V) larger than 0.05 mag, is, according to the framework 
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developped in this paper, due to an absence of scattered light contamination, 
because of the proximity of the star and of the obscuring material (section 7) . 

To conclude this section on the reddening of planetary nebulae, it is worth 
mentionning the article of Stasihska ct al. (1992) (and references therein), 
which compare the interstellar extinction of P.N. by means of photometry 
and a standard extinction law, and by means of radio measurements. There 
seems to be a systematic difference (always in the same way) when the red- 
denings found by the two methods are compared. This difference increases 
with reddening. According to Stasihska et al. (1992), it can only be due to 
a modification of the interstellar extinction law for clouds at larger distances 
than 2.2 kpc from the sun. At these distances, Ry = Ay/ E[B — V) would be 
smaller than in the vincinity of the sun. 

Strictly speaking the results of Stasihska et al. (1992) state that: 

• The observed extinction law depends on the distances {Iq and do) from the 
observer. 

• Standard extinction laws do not properly describe the interstellar reddening. 

Both items agree with the ideas discussed so far. That interstellar extinction 
varies with distance from the observer can hardly be explained by a modifica- 
tion of its properties, which would lead to a questioning case of geocentricism: 
why should the properties of interstellar grains vary with distance from the 
sun? 
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Table 3 

Sitko et al. (1981) shell stars 



star 




s - y (2) 


^(3) 






60 


(6) 


inn (7\ 
100* y ) 

lOOci 


BD+40° 4124 


B2 


0.77 


1.00 


b=+2.77 


108 


1.7 




0.9 


BD+61°154 


B8 


0.55 


0.66 


b=-0.95 


300 


5.0 




0.8 


V380 Ori 


AO 


0.48 


0.49 


4.528 


269 


4.5 




0.7 


HD259431 


B5 


0.28 


0.43 


b=+0.67 


290 


4.5 




0.5 


HD45677 


B2 


0.04 


0.28 


0.307 


355 


5.0 




0.2 


AB Aur 


B9-A0 


0.12 


0.15 


0.928 


144 


4.0 




0.3 


HD50138 


B8-B9 


+0.01 


0.10 


b=-3.14 


290 


4.5 




0.4 


HD190073 


AO 


0.09 


0.10 


0.118 


5000 


4.0 




0.7 


HD163296 


AlV 


0.07 


0.05 


b=+1.49 


122 


4.5 




0.9 



(1) Spectral type. 

(2) B — V from Sitko et al. (1981); updated with Simbad and the Lausanne Uni- 
versity database (http://obswww.unige.ch/gcpd). 

(3) E{B — V) deduced from spectral type and Fitzgerald (1970). The uncertainty 
on E{B — V) should not exceed ±0.05 mag.. 

(4) Large scale value of E{B — V) from Schlegel et al. (1998). For stars in the 
galactic plane, galactic latitude is given. 

(5) Estimated distance (pc) from the sun (Hipparcos). 

(6) Infrared color, around the star, of the interstellar cirrus. This is the 
slope of the linear regression between the 100 /xm and the 60 /xm images, on a 
2° X 2° field centered on the star. 

(7) Ratio of the 100 /xm emission at the position of the star, to what should be 
the 100 /im emission of the cirrus, estimated from the 60 /xm emission at the 
position of the star, and the linear regression found between the 100 /xm and 
60 /tm emission (item (6)). If the reddening of a star is mainly due to circum- 
stellar dust, should be small. On the contrary, if the extinction is due to 
foreground dust will be close to 1. 
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10 Stairs with circumstelleir dust; Sitko et al. (1981) shell stairs 

The first columns of Table 3 fist the main characteristics of the shell stars used 
in Sitko et al. (1981). In this paper, the authors show that when intrinsic and 

foreground reddenings can be separated, the bump feature associated with the 
intrinsic reddening is less than expected for a 'normal' star extinguished by 
the same amount of foreground interstellar matter. 

Stars in Table 3 are ordered by decreasing E(B — V). HD44179, the Red 
Rectangle, was excluded from the table, since we know today that we do not 
observe the star, hidden behind a thick torus of dust, but scattered starlight 
escaping from the edges of the torus (Waelkens et al., 1995). ZCMa (HD53179) 
and HD31648 were also discarded because the noise level of their LWP and/or 
SWP spectra is nearly as high as the signal itself. 

All stars in Table 3 have cirrus on their line of sight. Column Ed is the average 
opacity of the cirrus in the direction of each star (Schlegel et al., 1998) and 
column Jioo/Zeo its mean infrared color. Since the extinction due to these 
cirrus can be important, a major problem is to separate the reddening by 
circumstellar dust from the one due to the cirrus. 

10.1 Separation of intrinsic and foreground reddenings 

Prom the mean linear relation which exists between the 100 /xm and the 60 /xm 
emissions (seventh column of Table 3) in the vicinity of each star, we can cal- 
culate what should be the 100 fim emission of the cirrus at the position of 
the star, /loo.ci, for the 60 /im surface brightness measured by IRAS at the 
same position. The comparison of /ioo,ci to /loo,*, the measured 100 /xm sur- 
face brightness at the position of the star, gives an indication on the respec- 
tive proportions of foreground and intrinsic reddenings: if the proportion of 
foreground reddening is large /ioo,*/Aoo,cj should be close to 1, whfie a large 
proportion of intrinsic reddening should imply a high 60 fim surface brightness 
at the position of the star (because of the heating of the dust surrounding the 
star); thus a low /ioo,*/Aoo,a- 

Additional indications on the origin of the reddening of the stars will come 
from the distance; from the comparison of with the reddening of the stars; 
and from the level of this reddening, since intrinsic reddenings, as for P.N., 
can be presumed not to exceed 0.1-0.2 mag. 

These criteria applied to Table 3 permit the separation of BD+40°4124, BD+61°4124, 
V380 Ori, and HD190073, on the one hand, from HD45677, AB Aur, and 
HD50138 on the other. The cases of HD259431 and of HD163296, both at 
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low galactic latitude, arc more difficult to fix. HD259431 has an intermediate 
-^ioo,*/-^ioo,ci value, and its reddening is larger than what can be expected for 
a normal intrinsic one, which may be due to a better equilibrium between 
intrinsic and foreground extinctions. For HD163296 /ioo,*/-^ioo,ci is close to 1 
but this pre-main sequence star is close to the sun, and its reddening is low 
and may be due to the dust surrounding the star. 

The main source of reddening for the former stars, BD+40°4124, BD+61°4124, 
V380 Ori, and HD190073, is probably foreground dust. If so, the low fore- 
ground reddenings (<< 0.1 mag.) Sitko ct al. (1981) find for V380 Ori and for 
HD190073 are underestimated. This should also be the case for BD+61°4124, 
for which Sitko et al. (1981) give Ed ~ 0.3 mag.. The proportion of inrinsic 
reddening for HD45677, AB Aur, and HD50138, is larger. 



10.2 The bump feature and intrinsic reddening 



Sitko et al. (1981) tried to correct the extinction curves of AB Aur, HD50138, 
HD259431, and HD45677, for their foreground reddening, using a standard 
extinction law. This correction being done [Figm^c 9 of Sitko ct al. (1981)] the 
authors find a bump feature (attributed to intrinsic reddening) much lower 
than expected from a 'normal' reddening by foreground dust. 

This conclusion is supported by a direct analysis of the reduced spectrum of 
the stars (Figures 3 and A. 2, Table A.l). There are five stars in Table A.l 
with reddening less than 0.2 mag. The four stars whose reduced spectrum 
can be fitted without a scattered light component are all Sitko et al. (1981) 
shell stars. The bump-like depression in the spectrum of these stars, centered 
around 2400 A, is due to the Fe lines. HD50138 and AB Aur in particular have 
reddenings well above the 0.05 mag threshold. Their reduced spectrum can be 
compared to that of HD23480 (Figure 3), of lesser reddening, but with a clear 
2200 A bump feature. 

For HD259431 and HD45677, which have higher reddening, intrinsic and fore- 
ground extinctions are mixed and conclusions from the size of the bump are not 
so easy to reach. Comparison of their spectrum to HD62542 or to HD149757 
(C Oph) is meaningless since these two stars are surrounded by dust (Cardelli 
& Savage, 1988). 
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11 Effect of scattered steirlight contamination on distance esti- 
mates 

A component of scattered light in the spectrum of reddened stars, if it is not 
identified and separated from the direct starhght, will give an impression of 
less extinction by large grains than it really is the case. A star for which the 
extinction is underestimated appears farther, if its distance is estimated by 
means of photometric methods, than it really is. Consequently, photometric 
distances to reddened stars will be overestimated when compared to trigono- 
metric distances if the presence of scattered light is large enough to modify 
the estimate of the reddening of the star. 

This result can be retrieved through analytical considerations, since the pres- 
ence of scattered hght introduces a term which grows with 1/A: scattered light 
in the spectrum of reddened stars will always modify extinction curves by de- 
creasing the visible slope [oc E{B — V)], and therefore lead to an underestimate 
of the reddening. 

The difference between the reddening estimated from observation, E(B—V)obs: 
and the true E{B — V), deduced from the analysis of the spectrum of HD46223 
(Zagury, 2001a), is ~ 0.1 mag.. For the stars with high reddening used in 
Zagury (2001d) (Figure 4 of Zagury (2001d)), E{B - V) - E{B - V)obs can 
be up to 0.2 mag. The corresponding errors on the distance modulus [DM ~ 
5 log d — 5 ~ y — My — 3E(B — V),V and My are the observed and absolute 
magnitudes of a star in the K-band] are +0.3 mag and -1-0.6 mag, meaning 
that estimates of distance by means of UBV photometry can be out by up 
to Ad/d = +12%. This relative error can be enhanced if the reddening is 
estimated from wavelengths shorter than the V-band (H/9 for example), since 
scattered light contamination increases with decreasing wavelength. 

Hipparcos distances have been compared to photometric distances for some 
star clusters (see van Leeuwen & Hansen (1997); van Leeuwen (1999); Robi- 
chon et al. (1999); Stello & Nissen. (2001) and an article of Floor van Leeuwen 
at Hipparcos website), and for some planetary nebulae (Acker et al., 1998; Pot- 
tasch & Acker. 1998). There is agreement between the two methods for some 
objects, but there are cases for which the two methods lead to different results, 
with a systematic trend for Hipparcos distances to be smaller, by up to 10%, 
than the classical photometric methods. 

Since the same particularity exists for star clusters and for P.N., the cause is 
probably not linked to a specific class of objects (such as young clusters, as 
argued by van Leeuwen & Hansen (1997); see Stello & Nissen. (2001) who also 
rejects this hypothesis). It is not likely either that the reduction procedure 
of Hipparcos data be the cause of the discrepency; it would then affect all 
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Hipparcos data, while there is a good agreement between Hipparcos and other 
ground-based parallaxes (Harris et al., 1997). 

I wondered if the presence of scattered hght in the spectrum of reddened ob- 
jects could influence the determination of their reddening and distance, and 

contribute to the differences observed between Hipparcos and photometric dis- 
tances. In general, since extinction is much enhanced in the UV, UV extinction 
curves will be a powerful tool in determining the true extinction of a star for 
which the distance is to be determined by means of photometry. 



12 TowEirds a fit of the extinction curve 

12.1 The standard theory framework 

In the framework of the standard theory of extinction (i.e. if we do observe 
the direct light from reddened stars) the comparison of extinction curves, nor- 
malized to the same column density of large grains, is necessary to understand 
the relationships between the different types of particles responsible for the 
overall extinction. If the particles responsible for the UV extinction were the 
same large grains responsible for the visible extinction -or are in a constant 
proportion to these large grains- UV normalized extinction curves should su- 
perimpose well. In the contrary extinction curves normalized in the visible will 
differ in the UV. 

Although logical, this reasoning leads to a dead-end. Extinction curves nor- 
malized by E[B — V) have aleatory behaviors in the UV, meaning that the 
relations which clearly exist between the different parts of a non-normalized ex- 
tinction curve [Savage et al. (1985), or Zagury (2001d) for the relation between 
visible and far-UV extinctions] are lost when extinction curves are normalized 
by E{B — V)\ the number of parameters necessary to fit normalized extinction 
curves grows to 6, in contradiction with the empirical finding of Cardelli et 
al. (1989): 'The most important result presented here is that the entire mean 
eoctinction law, from the near-IR through the optical and lUE-accessihle UV, 
can he well represented by a mean relationship which depends upon a single 
parameter. . . The deviations of the observations from the mean relation (e.g.. 
Figs. 1 and 2) are impressively small.' [Cardelli et al. (1989), p. 252-253. The 
'mean relation' of their Fig. 1 hnearly relates, for UV wavelengths. Ax/ Ay and 
Ry^ = E{B - V)/Av, thus to E{B - V)]. 

Moreover the Fitzpatrick & Massa (1990) fits correctly reproduce the obser- 
vations in the UV but diverge from it in the visible. For three of the stars of 
Figure 3, the fit coefficients were calculated by Fitzpatrick & Massa (1988) 
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Fig. 3. Recapitulative plot of all the reduced spectra fitted in the Appendix. Fits are 
represented as dashes and reduced spectra as solid lines. The spectra are scaled by 
an arbitrary factor. E'ut,, is the reddening found from the fit, which can be slightly 
different from E{B — V) deduced from photometry (sec Table A.l). Fits are given 
in Table A.l. Sitko et al. (1981) stars are labeled with an 'S'. 
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Fig. 4. Fitzpatrick k. Massa (1990) fits (dashes) for three reduced spectra of Fig- 
ure A.l (also represented on Figure 3). The fits are good for the UV part of the 
spectrum (less in the bump region) but diverge in the visible. 
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or by Aiello et al. (1988). Figure 4 plots the reduced spectra of these stars 
along with their Fitzpatrick & Massa (1990) fit. The fits are good at UV wave- 
lengths (except in the bump region) but clearly diverge from observation in 
the visible. This absence of constraint of the Fitzpatrick & Massa (1990) fit by 
the visible part of the extinction curve is a first explanation to the absence of 
relationship found between the coefficients of the fit and E(B — V) (Jenniskens 
& Greenberg, 1993; Aiello et al., 1988). 



In a recent paper (Fitzpatrick, 2004), E. Fitzpatrick recognizes the limitations 
of the Fitzpatrick & Massa (1990) fit, and argues that the correlations found 
by Cardelli et al. (1989) between the different parts of non-normahzed extinc- 
tion curves is only due to a too small number of stars [29 stars are used by 
Cardelli et al. (1989)]. This argument doesn't stand. Firstly, Cardelli et al. 
(1989) stars are field stars, which were taken from Fitzpatrick & Massa (1988) 
and had available near-IR photometry at the time the paper was written, with 
no particular assumption governing their choice. Why shouldn't the following 
fifty or hundred stars follow the same relationships, except for possible de- 
viations due to uncertainties on spectral types? Secondly, even the 29 stars 
of the Cardelli et al. (1989) study present strong divergences in the UV if 
they are normalized by E{B — V) (see Figure 6 of Bless & Savage (1972)). 
Lastly, the existence of a relationship, between the visible and far-UV parts of 
extinction curves, which is evident in the study of Cardelli et al. (1989), was 
independently retrieved by Zagury (2001d), on the basis of the physical ideas 
developed in the present paper. 



It is therefore the normalization process of extinction curves which is to be 
questioned, and not the number of stars used by Cardelli et al. (1989). 



If scattered starlight is added to direct starlight, the light we receive from a 
reddened star will be composed of direct starlight, extinguished by large grains 
in proportion of the column density (thus proportional to E{B — V)), on the 
one hand, and of the starlight scattered by hydrogen, also extinguished by large 
grains, but proportional to the square of the column density [E{B — V^], on 
the other. Because of this difference in the dependence of the two components 
of the extinction curve on E{B—V) , normalization by E{B—V) is meaningless. 
It can only increase the differences between extinction curves, and mask their 
relationships. The normalization by E{B — V) is the second reason, in addition 
to its limited wavelength domain of application, why no relation was ever 
found between the far-UV and visible parts of the Fitzpatrick & Massa (1990) 
decomposition of the extinction curve. 
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Fig. 5. Reddening estimated from the UV reduced spectra, Euv, versus the reddening 
estimated from the visible photometry, Eyis, for Sitko et al. (1981) (o) and Zagury 
(2000b) stars (*). 

12.2 A first step towards the fit of the eoctinction curve 



Rayleigh scattering by hydrogen, with an additional extinction by large grains, 
should be fitted by a function proportional to (Rayleigh, 1871): 



1 _ a ri-E.s 



(17) 



The e^"^^"/^ term is for the extinction by large grains of the scattered light. 
The e^^/'^ term accounts for extinction by the gas. g should grow with the 
gas column density. This term is usually considered as small compared with 
the extinction by large grains, but it may gain in importance in the UV. It 
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Fig. 6. PE^^, the coefficient which defines the proportion of scattered hght, versus 
Euv (equation 19) for Sitko et al. (1981) (*) and Zagury (2000b) (o) stars. PE^^ 
roughly behaves as E^^. 

must affect direct starlight as well as scattered light. 



Outside the bump region, the fit of the reduced spectrum of a reddened star 
should ressemble: 



-.-2^-% , t:, 2 r/^\-2^ --2- 



e-^^e-^ + E/f{j)e-'^e-t^, (18) 



A in //m. 

Since direct and scattered light cross media which should be very much alike, 
we can suppose Qd = 9 and = Eg = E^v {E^v is the true reddening, 
determined from the UV part of the reduced spectrum, which should be close 
to E{B — V)obs deduced from photometry and spectral type). 
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Fig. 7. g (equation 19) versus E^v for Sitko et al. (1981) (*) and Zagury (2000b) 
(o) stars. The fits are probably not precise enough to fix the exact dependence of g 
on Euv [or E{B - V)]. The straight hne gives g ~ 0.00055£;„^. 

/(1/A) contains the dependence of the scattered light, corrected for the ex- 
tinction by large grains, on 1/A. / also depends on the distances to the cloud 
and to the star. 

Observations I have used (Zagury, 2001a, d), give an / dependence on A as 
1/A^. Therefore, the analytical expression of interstellar extinction outside 
the bump region should look like: 

e-^^e-*(l + /?^) (19) 
P depends on the distances to the cloud and to the star. 
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Fig. 8. The UV flux of the unreddened AO star HD104878 is compared to the flux 
of AB Aur multiplied by e^'^^^^ (which accounts for the reddening of AB Aur) and 
to the flux of HD190073 multiplied by e°-25/^. The spectra superimpose well in the 
visible (not represented here, see Figure A. 2 for AB Aur) and in the near-UV. In 
the far UV the spectra are difficult to compare because of the large Fe and Lya 
absorptions. The positions of the Fe line and of the 2200 A bump are marked by 
vertical lines. 

Figure 3 summarises all the fits calculated (from expression 19) and included 
in the appendix. The spectra are ordered from top to bottom by decreasing 
reddening. The expression of each fit is in Table A.l. Although there are 
uncertainties on these fits and they probably can be improved, the coefficients 
found already give an idea of the relationships which may exist between the 
parameters of expression 19. 

Figure 5 confirms the relationship, already evidenced in Zagury (2001d), be- 
tween the visible and the far-UV parts of the extinction curve. E'^^, the value 
of the reddening in the direction of each star obtained from the analysis of the 
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extinction curve, is plotted against Ems, the reddening determined from the 
star photometry. There is good agreement between the two values. 

The PEuv'^ coefficient which fixes the importance of the scattered light is 
plotted against E^v in Figure 6. It roughly behaves as E^y'^, as expected, with 
P ~ 0.035 (neglecting the dependence on the geometry). 

The gas extinction term g affects the far-UV part of the fits, permitting a 
steeper curve in the far-UV. It is not clear whether the steepness of the far-UV 
extinction curve is due to this term or to the wing of the Lya absorption line. 
AB Aur (Figure A. 2) and HD190073 for instance have steep far-UV reduced 
spectra. Since the reddening in these two directions is low it can be surmised 
that gas extinction is negligible -even in the UV- so that the steepness of 
the far-UV extinction curve in these directions is due to the wing of the Lya 
absorption. Comparison with the non-rcddenned star HD104878 shows this is 
certainly the case for AB Aur, and probably also for HD 190073 (Figure 8). 
The absence of a clear relation in Figure 7 plot of g versus Euv will also be 
related to this difficulty of a precise determination of g. The pertinence of a 
gas extinction term can probably also be investigated from a theoretical point 
of view by the comparison of the extinction cross-section of the gas to the 
values of g {a few 10~^, if A is in /im) given in Table A.l. 



13 Conclusion 

The hypothesis of the standard theory of extinction, that the UV fight we 
receive from a reddened star is mainly the direct light from the star, and its 
corollary, the three component interstellar grain model, are not borne out by 
observation. 

The only alternative to the standard theory is that the spectrum of a reddened 
star is contaminated by light scattered at small angular distances from the 
star. The far-UV part of the extinction curve is observed to be less than the 
continuation of the visible extinction by large grains because of the importance 
of the scattered light component in the UV. 

There is no need for particles such as PAH to explain the far-UV light we 
receive from a reddened star, with the probable implication that these particles 
do not exist in the interstellar medium. This is not so much surprising since 
one can wonder how particles we are not even able to synthesize on earth, can 
proliferate in the cold and low density environment of an average interstellar 
cloud. 

In the first part of the paper I have discussed several aspects of coherent scat- 
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tering in the near forward directions, which is the unique reason I have found 
to explain the importance of the scattered hght component in the spectrum 
of reddened stars (Zagury, 2001b). There is up to date no theoretical treat- 
ment of this phenomenom, only briefly mentioned in some textbooks as the 
Bohren & Huffman (1983). The reason might be the specific conditions nec- 
essary to its observation, which require large distances between the source of 
light, the scattering medium, and the observer; and a non-homogenous (tur- 
bulent) medium. 

Elementary calculations show that coherent forward scattering by hydrogen 
easily accounts for the order of magnitude found for the scattered light com- 
ponent of the spectrum of reddened stars, and will give the observed 1/A^ 
wavelength dependence. This scattering occurs within angles of 10^*^" or less 
from the direction of the stars, which represents a surface of radius a few 
100 km for a cloud 100 pc away. It imphes that interstellar clouds still have 
structure at scales of a few hundred km. 

If the cloud is close to the star {do = 0) variations of the pathlengths for 
the scattered light will be large (compared to UV wavelengths) and random, 
even at very close angular distance from the star. The effect of coherence will 
be annihilated. This latter result was announced by Sitko et al. (1981), who 
showed that hot stars with circumstellar dust have a 2200 A bump smaller 
than expected from their reddening. Sitko et al. (1981) result was extended 
to planetary nebulae (section 9). It is concluded that circumstellar matter 
participates httle or nothing to the 2200 A bump. The bump, when present, 
is mainly due to a foreground cloud. The extinction curves for the Sitko et al. 
(1981) stars with 0.1 < E{B — V)< 0.2 are a straight line from the visible to 
the UV: we observe the extinction of the direct starlight by the large grains 
(with a continuous extinction law) only. 

The diminution of the 2200 A bump when extinction comes from circumstellar 
material emphasizes the link between the coherent scattered light and the 
bump. It confirms an opinion I have already expressed (Zagury, 2001b, 2002b) 
that the bump is not an absorption feature but an interruption of the scattered 
light, or, as I believe today, a compensation between (1) the extinction by the 
gas of the direct starlight, and (2) the scattered hght. 

Standard extinction laws tend to underestimate the true extinction by large 
grains. If the distance of a star is estimated by means of photometry, the 
unreddened flux of the star will be underestimated, and the star will be be- 
lieved to be farther off than it really is. The mismatch can be larger than 10 %, 
and should be observed when photometric distances of reddened stars are com- 
pared to parallaxes. This may be related to on-going problems (van Leeuwen, 
2004) on the comparison between Hipparcos results to distances determined 
by photometric means. 
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Based upon the dependence on A of the far-UV extinction curve in a few 
directions found in Zagury (2001a,d), I have discussed what should be the 
fit of the extinction curve (outside the bump region). This early attempt to 
give a physical meaning to the fit of the extinction curve was applied to the 
stars used in this paper and in Zagury (2000b). The relationship between far- 
UV and visible extinctions found by Zagury (2001d) is well confirmed. The 
coefficient which determines the importance of the scattered light component 
roughly behaves as Nfj, as expected. The fit used in section 12 and in the 
appendix introduces a gas extinction term oc e~^/^ which accounts in part 
for the steepness of the extinction in the far-UV. It is not clear whether this 
effect is real or due to the wings of the Lya absorption. 

In section 12.1 1 have also discussed the reasons why no relation was ever found 
using the Fitzpatrick & Massa (1990) fit of the extinction curve, and more 
generally why this fit is not appropriate -despite a new attempt, announced in 
Fitzpatrick (2004), to model simultaneously the spectra of the reddened and 
standard stars- to a correct representation of the extinction curve. 

A last imphcation of this work, foreseen, and refused, by Fitzpatrick (2004), is 
that the true interstellar extinction law, and the average extinction properties 
of interstellar grains, are the same in all directions. 



A Appendix: Fit of extinction curves 

The reduced spectra of the stars used in Zagmy (2000b) are fitted in Fig- 
ure A.l, using the general formulation of section 12, equation 19. There is not 
absolute uniqueness of each fit, and the relevance of the gas extinction param- 
eter g still remains to be proved (section 12). In absence of more constraints I 
tried to determine a good fit to each of the spectra. The fits extend from the 
visible to the far-UV, the visible spectra are replaced by the photometric UBV 
points retrieved from the Lausanne database (http:/ /obswww.unige.ch/gcpd). 

The reduced spectra of Sitko et al. (1981) stars can be classified, as I have 
done in previous papers, according to how far towards the visible the scattered 
component can be perceived. HD163296, HD190073, HD50138, AB Aur, follow 
the same extinction law in the visible and in the UV (Figure A. 2), thus have a 
negligible scattered light contamination. For HD45677, of moderate reddening, 
the visible extinction law extends down to the bump region, the scattered 
component appears as superimposed on the tail of the far-UV extinction by 
large grains. The reduced spectrum of each of the four stars of higher reddening 
(HD259431, V380 Ori, BD+61° 154, BD+40°4124) muhiplied by A^ is an 
exponential of 1/A in the far-UV (Figure A. 3); the exponent, 2Euy, is close to 
2E{B - V) (Figure 5). The far-UV slope of BD+40°4124, the star with the 
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2 4 6 8 10 

Fig. A.l. Fit (dashes), outside the bump region, of the reduced spectra (sohd hues) 
of Zagury (2000b), with an arbitrary logarithmic y-scahng. The visible part of the 
spectra is represented by the phbotometric UBV points. Dotted lines arc the direct 
and scattered starlight components. The expression of each fit can be found in 
Table A.l 
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highest reddening, is difficult to determine, because of the small wavenumber 
interval between the bump on the short wavenumber side, and the Lja line 
on the other side. Figure A. 2 plots, as it was done for the stars of Zagury 
(2000b), the reduced spectra of Sitko et al. (1981) stars (except for HD190073, 
see figure 8) and a possible fit deduced from equation 19. 
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Fig. A. 2. Fits, outside the bump region, for the reduced spectra of Sitko et al. 
(1981) stars. Fit are represented as dashes and reduced spectra as sohd hues, y-axis 
is logarithmic. The spectra are scaled by an arbitrary factor. E^v is the reddening 
found from the fit, which can be slightly different from E(B — V) deduced from 
photometry (see Table A.l). The analytical expression of each fit is in Table A.l. 
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Fig. A. 3. Far-UV reduced spectra of Sitko et al. (1981) stars of high reddening, 
multiphed by A^. Each spectrum is well fitted by an exponential of 1/A, with an 
exponent 2Euv close to 2E{B - V) (Figure 5). Only for BD+40° 4124 is the slope 
difficult to determine because of the importance of the 2200 A bump and of the Lya 
absorption. 
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Table A.l 

Fits used in Figures 3, A.l, and A. 2 



star 



E (3) 



fit 



BD+61°154 

HD200775 

V380 Ori 

HD62542 

HD259431 

HD45677 

HD37903 

HD149757 

AB Aur 

HD50138 

HD190073(6) 

HD23480 

HD163296 



HD23753 
HD58050 
HD104878 
HD32630 
HD199081 
HD31726 
HD74273 
HD214680 
HD104878 
HD23753 
HD 104878 
HD215573 
HD119765 



0.66 
0.55 
0.49 
0.38 
0.43 
0.28 
0.32 
0.33 
0.15 
0.10 
0.10 
0.09 
0.05 



0.62 
0.60 
0.47 
0.43 
0.37 
0.33 
0.33 
0.33 
0.15 
0.14 
0.12 
0.10 
0.09 



b=-0.95 
11.946 
4.528 
0.76 
b=+0.67 
0.31 
2.02 
0.60 
0.928 
b=-3.14 
0.118 
0.64 
b=+1.49 



1.25 . 

oc e A e 



0.00032 / „ 

fl + Mi 



1.2 _ 0.0002 

OC e A e a4 



0.94 0.0(1(12 
3 



OC e A e A 

0.86 0.0005 

OC e A e a4 

0.74 0.00019 



OC e A e 



A* 



0.66 _ 0.00035 

OC e A e a4 

0.65 _ 0.00022 

OC e A e A* 

0.66 0.00021 

OC e A e A* 



l + Mll) 
1 , 0.0031 A 



1 + 



0.0067 \ 

1 I aoo47A 
^ + ^?-) 

1 I a0045\ 

1 I aoo2oA 



0.31 

OC e A 

0.28 

OC e A 

0.25 

OC e A 



20 _ 0.00016 

OC e A e a4 



^-|^ _l_ 0.00046 \ 



0.17 

OC e A 



S 
uv2 

S 
uv2 

S 

S 
uv2 
uv2 

S 

S 

s 

uv2 
S 



(1) Standard star used to establish the reduced spectrum of the star. Each of these 
reference star was corrected for its slight reddening (see Table A. 2) according 
to Zagury (2001c)). 

(2) E{B — V) from UBV photometry and spectral type. Corresponds to Eyig in 
Figure 5. 

(3) Euv deduced from the UV part of the reduced spectrum. 

(4) Large scale value of E{B — V) from Schlegel et al. (1998). Galactic latitude 
for stars close to the galactic plane. 

(5) 'S' is for Sitko et al. (1981) sheh stars, uv2 for Zagury (2000b). 

(6) HD190073 was not represented on Figures A.2 and 3, but on Figure 8. 
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Table A.2 
Standard stars 



star 


S.T.W 


5 - y (2) 


^(3) 


HD 104878 


AOV 


-0.011 


0.00 


HD119765 


AlV 


+0.007 


-0.01 


HD199081 


B5V 


-0.136 


0.02 


HD2 14680 


09V 


-0.205 


0.10 


HD215573 


B6IV 


-0.123 


0.02 


HD23753 


B8V 


-0.070 


0.04 


HD31726 


B2V 


-0.212 


0.03 


HD32630 


B3V 


-0.146 


0.05 


HD58050 


B2Ve 


-0.187 


0.05 


HD74273 


B1.5V 


-0.209 


0.04 



(1) Spectral type. 

(2) B—V from the Lausanne University database (littp://obswww.unige.cli/gcpd). 

(3) E{B — V) deduced from spectral type and Fitzgerald (1970). 
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